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ABSTRACT 

In this paper we present a new method to extract cosmological parameters using 
the radial scale of the Baryon Acoustic Oscillations as a standard ruler in deep galaxy 
surveys. The method consists in an empirical parametrization of the radial 2-point cor- 
relation function, which provides a robust and precise extraction of the sound horizon 
scale. Moreover, it uses data from galaxy surveys in a manner that is fully cosmology 
independent and therefore, unbiased. A study of the main systematic errors and the 
validation of the method in cosmological simulations are also presented, showing that 
the measurement is limited only by cosmic variance. We then study the full informa- 
tion contained in the Baryon Acoustic Oscillations, obtaining that the combination of 
the radial and angular determinations of this scale is a very sensitive probe of cosmo- 
logical parameters, able to set strong constraints on the dark energy properties, even 
without combining it with any other probe. 

Key virords: data analysis - cosmological parameters ~ dark energy - large-scale 
structure of the universe 



1 INTRODUCTION 

Finding the physical origin of the accelerated expansion of 
the Universe is one of the most important scientific problems 
of our time, and is driving important advances in XXIst cen- 
tury cosmology. Several observational probes to study the 
nature of the mysterious dark energy, which powers that 
expansion, have been proposed. Among them, the measure- 
ment of the scale of the Baryon Acoustic Oscillations (BAO) 
in the galaxy power spectrum as a function of redshift is 
one of the most robust, since it is insensitive to system- 
atic uncertainties related to the astrophysical properties of 
the galaxies. Moreover, it provides information about dark 
energy from two different sources: the angular diameter dis- 
tance, through the measurement of the BAO scale in the 
angular distribution of galaxies, and the Hubble parameter, 
through the measurement of the BAO scale in the radial 
distribution of galaxies. 

There are some measurements of the BAO scale in the 
purely radial direction ( ,Gaztanaga et al.^^2009( [Kazin et al.| 



2010 1, but most of them use mainly the monopole of the 3-D 



correlation function ( 


Eisenstein et al.||2005 


Hiitsi||2006a|b| 


Percival et al.||2007t 


Padmanabhan et al. 


20071 Okumural 


et al.|2008 


Sanchez et al.|2009 


Anderson et al.|2012 


1. This 



is not the most sensitive determination of the BAO scale, 
and in fact is strongly model-dependent. 



In this paper we propose a new method to extract the 
evolution of the radial BAO scale with the redshift, and 
explain how to use it as a standard ruler to determine cos- 
mological parameters. We use data from galaxy surveys in 
a manner that is fully cosmology independent, since only 
observable quantities are used in the analysis and therefore, 
results are unbiased. A method based on the same idea for 
the measurement of the angular BAO scale was described 
in 



Sanchez et al. (20111 and provided the measurement of 



(20121 
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the angular BAO scale presented in Carnero et al 
Here we present how to extract the radial BAO scale using 
a similar approach. 



The method is designed to be used as a strict standard 
ruler, and provides the radial BAO scale as a function of the 
redshift, but we do not intend to give a full description of 
the radial correlation function. This approach is more robust 
against systematic effects, and in fact we demonstrate that 
the measurement is only limited by cosmic variance, since 
the associated systematic errors are much smaller than the 
purely statistical errors. 
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2 GALAXY CLUSTERING AND 
OBSERVABLES 

One of the main statistical probes of the properties of the 
matter distribution in the Universe is the 2-point correlation 
function, ^(r), which is defined as the excess joint probabil- 
ity that two point sources {e.g. galaxies) are found in two 
volume elements dVidV2 separated by a distance r compared 
to a homogeneous Poisson distribution ( Peebles|1980 1. If the 
fluctuations on the matter density field are Gaussian, this 
function contains all the information about the large scale 
structure of the Universe. 

However, what is observationally accesible is the distri- 
bution of galaxies in angle-redshift space, not directly the 
matter distribution in real space. For each galaxy we de- 
termine its angular position in the sky and its redshift. To 
obtain ^ (r) we need to convert the measured redshift to a co- 
moving distance, for which a cosmological model is needed. 
Therefore, the 3-D correlation function is not observable in a 
cosmology independent way for a given galaxy survey. More- 
over, the observational techniques to obtain the angular po- 
sition in the sky and the redshift are completely different 
and independent. Consequently, the most natural way of 
measuring BAO from the observational viewpoint is, on the 
one hand, the measurement of the angular correlation func- 
tion as a function of the angular separation of galaxy pairs, 
and on the other hand, the radial correlation function as a 
function of the redshift separation of galaxy pairs, and then 
extract the BAO scale from each function. Both of them are 
defined in complete analogy with the 3-D function, as the 
excess of probability with respect to a homogeneous Poisson 
distribution, and both are observable in a cosmology inde- 
pendent way, providing, together, information about the full 
3-D large scale structure that can be described in any cosmo- 
logical model. These aspects have been recently pointed out 
in [Bonvin &: DurreT] ( |2011| | , and here we present a practical 
implementation of a cosmology independent measurement, 
in this case for the radial BAO scale. 

Several effects must be taken into account when the 
large scale structure of the Universe is studied from galaxy 
surveys, like bias, redshift space distortions and non-linear 
corrections. In this paper we take all of them into account, 
both in the theoretical calculations and in the observational 
analysis. For this particular measurement, the influence of 
the non-linear effects is relatively small, since we are analyz- 
ing large scales, around the BAO scale. Also the influence 
of bias is small, because the analysis method is designed 
to minimize its impact on the final result. However, redshift 
space distortions are very important, since we are measuring 
along the line of sight. 



3 THE SOUND HORIZON SCALE 

The BAO are a consequence of the competition between 
gravitational attraction and gas pressure in the primordial 
plasma, which produced sound waves. Differences in density 
created by these sound waves leave a relic signal in the sta- 
tistical distribution of matter in the Universe, defining a pre- 
ferred scale, the sound horizon scale, which is a robust stan- 
dard ruler from which the expansion history of the Universe 
can be inferred. The BAO scale is very large, ~ 100 Mpc/h, 



which poses an important challenge for observations, since 
surveys must cover large volumes to map such a distance. 
However, the very large scale of BAO also has certain advan- 
tages, because structure formation at these scales is rather 
well understood, and details about the description of galaxy 
formation and astrophysics do not compromise the accu- 
rate measurements of the standard ruler. In practice, the 
strength of the BAO standard ruler method relies on the 
potential to relate the position of the acoustic peak in the 
correlation function of galaxies to the sound horizon scale 
at decoupling. 

The current measurements of the BAO scale typically 
use the monopole contribution of the 3-D correlation func- 
tion. This is a mixture of the angular and radial scales, and 
therefore, does not contain all the information that can be 
extracted from the data samples. On top of that, all analy- 
ses have been done using what is usually called a "fiducial 
model" in order to determine distances. Rather than mea- 
surements, they should be understood as consistency checks. 
Considering this way of analyzing data, what the current re- 
sults have achieved is a verification of the consistency of the 
data with the standard ACDM model, but the analysis strat- 
egy is very limited if one wants to test other cosmological 
models against data, or at least, the whole analysis must 
be repeated from the very begining for every cosmological 
model one tries to test. 

On the other hand, if both the angular and the radial 
BAO scales are determined, which in principle is already 
possible with current data sets, the cosmological constraints 
that can be set using only the BAO scale are much stronger 
than the current values, since the radial and angular scales 
are sensitive to cosmological parameters in quite different 
ways, breaking degeneracies and allowing a very precise de- 
termination when combined. 



4 METHOD TO MEASURE THE RADIAL BAO 
SCALE 

The method we propose relies on a parametrization of the 
radial correlation function which allows extracting the ra- 
dial BAO scale with precision. It is based on a similar ap- 
proach previously proposed for the angular correlation func- 
tion (Sanchez et al. 20111. The suggested parametrization 



was inspired by the application of the Kaiser effect to the 
shape we used for the angular correlation function. 

It is important to remark that the calculations we per- 
form to check the goodness of the parametrization do in- 
clude redshift space distortions using the Kaiser description 



(IKaiser 19871 and also non-linearities using the RPT ap- 
proach ( Crocce & Scoccimarro '2006jl , and we expect the 
result to be correct for large scales, where the neglected ef- 
fects (fingers of god, mode-mode coupling) are very small. 
We have used a bias parameter 6=1 throughout the calcula- 
tion. The possible infiuence of bias is studied in the system- 
atic errors section. We use the linear power spectrum from 
CAMB ( Lewis et al.||2000 l, on top of which non-linearities 
and redshift space distortions are then included. 

The full recipe to obtain the radial BAO scale as a func- 
tion of redshift for a galaxy survey is as follows: 

(i) Divide the full galaxy sample in redshift bins. 

(ii) Divide each redshift bin in angular pixels. 
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Figure 1. Some examples of the description of the different cosmological models by the proposed parametrization. Models do include 
redshift space distortions and non-linearities using the RPT scheme. The parametrization is a very good description of all models at all 
redshifts, even if the considered errors are much smaller than the cosmic variance for the full sky. Errors are smaller than the dot size. 
We have used this level of precision to ensure that the systematic errors associated to the theory (non-linearities, bias) are small. 



(iii) Compute the radial correlation function by stacking 
galaxy pairs in each angular pixel, but do not mix galaxies 
in diflterent angular pixels. 

(iv) Parametrize the correlation function using the ex- 
pression: 

Cii (Az) =A + i?e-^^- - De-^"^^ + Fe- '^'"tf (i) 

and perform a fit to ^||(A2:) with free parameters A, B, C, 
D, E, F, Azbao and a. 

(v) The radial BAO scale is given by the parameter 
Azbao- The BAO scale as a function of the redshift is the 
only parameter needed to apply the standard ruler method. 
The cosmological interpretation of the other parameters is 
limited, since this is an empirical description, valid only in 
the neighborhood of the BAO peak. 

(vi) Fit cosmological parameters to the evolution of 
Azbao with z. 

In order to test the goodness of this parametrization, 
we have computed the radial correlation function for the 
14 cosmological models described in Table [l] in a redshift 
range from 0.2 to 1.5, always including redshift space dis- 
tortions and non-linearities as explained before. Then, we 
have applied the method to each model and each redshift. 
The parametrization describes the theory very accurately, 
since the values of the jndo] are close to 1, and the prob- 
abilities of the fit lie between 0.98 and 1.00 when the error in 
each point of the correlation function is arbitrarily fixed to 
2% (or ~ 10~^) for all bin widths and cosmological models. 
This error corresponds to a precision much better than the 
cosmic variance for the full sky, for all models in the full red- 
shift range, but the parametrization is able to recover the 
correct radial BAO scale for the 14 cosmologies. We have 
used this level of precision because the systematic errors 



h fijvj Hi, Qfe uiQ Wa ris 

0.70 025 0.044 OOO -1.00 0.95 
0.68 
0.72 

0.20 
0.30 

0.040 
0.048 

+0.01 
-0.01 

-0.90 
-1.10 

-0.1 
-1-0.1 

1.00 

Table 1. Summary of the 14 cosmological models used to test the 
method. Where empty, the fiducial values (first line) are assumed. 



coming from theoretical effects (non-linearities, bias, fingers 
of god) only affect this method if they induce a disagreement 
between the models and the parametrization, giving then a 
wrong measurement of the BAO scale. If the description is 
good to this level of precisionwe guarantee a small contribu- 
tion from these systematic effects, as will be shown in the 
systematic errors section below. The calculation of the real- 
istic errors is described in the next section. Some examples 
of these descriptions can be seen in Figure [l] 

Once we have verified the parametrization on the theo- 
retical predictions, we are ready to apply it in more realistic 
environments. For this purpose, in the next section we ap- 
ply this algorithm to a galaxy catalogue obtained from a 
N-body cosmological simulation. We will also compute the 
main systematic errors associated to this method. 
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Redshift 



Figure 2. N{z) in tiie used MICE catalogue. Tiie simulation 
contains 55 million galaxies in the redshift range 0.1 < z < 1.5. 
The vertical dashed lines show the limits of the redshift bins used 
in the analysis. 



5 APPLICATION TO A SIMULATED GALAXY 
SURVEY 

We have tested the method to recover the radial BAG scale 
using a large N-body simulation capable of reproducing the 
geometry (e.g. area, density and depth) and general features 
of a large galaxy survey. The simulated data were kindly pro- 
vided by the MICE project team, and consisted of a distri- 
bution of dark matter particles (galaxies, from now on) with 
the cosmological parameters fixed to the fiducial model of 
Table 111 The redshift distribution of the galaxies is shown in 
Figure |2| The simulation covers 1/8 of the full sky (around 
5000 square degrees) in the redshift range 0.1 < z < 1.5, and 
contains 55 million galaxies in the lightcone. This data was 
obtained from one of the largest N-body simulations com- 
pleted to date with comoving size Ltox ~ 3072h~^ Mpc 
and more than 8 x 10^ particles (nip = 2.3 x IO^/i'^Mq). 
More details about this simulation can be found in IFosalbal 
et al.l (|2008|) and iCrocce et al.l (|2010|). The simulated cat- 



alogue contains the effect of the redshift space distortions, 
fundamental for the study of the radial BAG scale. 

Data with similar characteristics will be obtained in fu- 
ture large spectroscopic surveys, such as DESpec ( |Abdalla| 



^ Laureijs et al.|20TT l. 



et al. 20121, BigBGSS (Schlegel et al. 20111 or EUCLID 



It is important to note that the radial BAG determina- 
tion needs a very large survey volume. We tried to extract 
the BAG peak from catalogs with smaller areas (200, 500 
and 1000 sq-deg), finding a very small significance (or no 
detection at all) in most cases. This is due to the fact that 
the statistical error related to the cosmic variance is spe- 
cially large for the radial correlation function, and therefore 
it can only be reduced by increasing the volume explored. 

Thus, we have divided the simulation into 4 redshift 
bins. We apply the method described in the previous sec- 
tion for each bin and obtain the correlation functions using 
the Landy-Szalay ( Landy Sz Szalay|1993 1 estimator. Fits are 



shown in Figure [3] We generate the random catalogues tak- 
ing into account the N{z) distribution in order to obtain 
the correct determination of the radial correlation function. 
This is a very important effect, since we are measuring in 
the radial direction and the distributions on the redshift 
coordinate have a large influence in measurements. The sta- 
tistical significance of the BAG observation in the flrst bin 
is very low (~ 1.4a) and it is consequently not considered 
in the cosmology analysis. Following the same approach of 
Sanchez et al. (20111, we have computed the statistical sig- 



nificance of the detection by measuring how different from 
zero the F parameter of the fit is (Eq.[T]), using its statistical 
error. 

We have computed a theoretical estimation of the co- 
variance matrix for the 3-D correlation function (Appendix 
[A| ), which is then applied to the estimation of the covariance 
for the radial correlation function, taking the corresponding 
value for the line of sight direction. This estimation relies 
on the model used for the simulation, but we expect a small 
variation of the error with the cosmological model. In any 
case, to keep the method fully model independent we have 
validated the calculation of the covariance matrix obtaining 
it by an alternative method directly from the simulation, 
which can also be applied to any real catalogue. We have 
used many realizations, dividing the total area in patches of 
different sizes and computing the corresponding dispersions 
and covariance matrices. We have then obtained the error 
for the total area, scaling this estimates to the full area of 
the measurement. Both determinations agree in the region of 
interest for the BAG scale measurement, as is shown in Fig- 
ure |4] There is a disagreement for small scales, which is com- 
ing from the incomplete description of the non-linearities in 
the theoretical calculation, where the mode-mode coupling 
effects are neglected and from boundary effects in the real- 
izations. 

It is important to remark that the volume which the 
radial BAG determination needs is very large. We are not 
able to obtain a significant observation for the first redshift 
bin and the contribution to the total error of the cosmic 
variance and the shot noise is comparable for all the other 
bins, which is shown in Figur e [4] A s already pointed out 
and verified by Crocce et al. (20111, in order to obtain a 



^ This catalogue is publicly available at http:/ /www. ice. cat/mice 



correct description of the error we need to sum linearly the 
cosmic variance and the shot noise contribution. This is not 
surprising, since galaxy surveys are designed to have a de- 
termined and fixed galaxy density, and this fact correlates 
both contributions to the error. 

The measured values of the BAG scale as a function of 
the redshift can be seen in Figure[5]as dots. The results for 3 
alternative bins shifted with respect to the nominal ones are 
also shown as triangles. These have not been used to obtain 
cosmology results, since they are, obviously, fully correlated 
with the red ones and are only shown for illustration and 
verification of the correct behaviour of the analysis method. 
We have used in this analysis the center of the redshift bin to 
obtain the prediction of the model, although what is really 
observed is the average within the bin. We have verified that 
they are very close if the N{z) distribution is smooth, as it 
is in this case. 
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Figure 3. Radial correlation functions measured in the MICE 
simulation for the 4 redshift bins described in the text (dots) 
compared with the proposed parametrization (solid line) . All the 
fits are good. The statistical significance of the BAO detection in 
the first bin is very low, and it is not used to set cosmological 
constraints. 
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Figure 4. Comparison of the different estimates of the error in 
the radial correlation function are shown. The estimates are MC 
samples (dots) and theoretical calculation (solid line) . They agree 
in the region of interest for the BAO analysis. The disagreement 
at low scales comes from the incomplete description of the non- 
linearities, where the mode-mode coupling effect is neglected, and 
from some boundary effects in the realizations, but does not affect 
the BAO scale measurement since it is outside the fitting region. 
The different contributions to the total error are shown. They 
come from Poisson shot noise (dotted line) and cosmic variance 
(dashed line). The contribution of the Poisson shot noise is not 
negligible. 




Figure 5. Measured radial BAO scale as a function of the redshift 
in the MICE simulation using the proposed method. Dots are the 
nominal bins and triangles correspond to displaced bins, and are 
measured only as a cross-check. All measurements are in good 
agreement with the theoretical prediction (solid line). 



5.1 Systematic Errors 

We have studied the main systematic errors that affect the 
determination of the radial BAO scale using this method. 
We have found a specific systematic effect associated to the 
algorithm, which is coming from the size of the angular pixel. 
Other systematic errors are generic and will be present in 
any determination of the BAO scale, namely, the infiuence 
of the non-linearities, the starting and end point of the fit 
to the correlation function and the possible influence of the 
galaxy bias in the measurement. 



6 E. Sanchez et al. 



0.35 




• 


1 sa-dea px 


0.3 




+ 


0.25 sq-deg px 


0.25 




■ 


0.0625 sq-deg px 


0.2 


t ^ 




0.0309 sq-deg px 


0.15 








0.1 








0.05 













. . . 1 


1 



10' 



10-^ 



Az 



Figure 6. Effect of the angular pixel size on the radial correlation 
function. The finite size of the angular pixel induces a change on 
the slope of the correaltion function at small scales. The effect 
is shown for pixels of sizes 0.0309 (stars), 0.0625 (squares), 0.25 
(crosses) and 1 (dots) degrees. The change in the slope arises from 
the smoothing effect produced by the inclusion in the calculation 
of galaxy pairs which are not exactly coUinear. This effect does 
not affect the determination of the BAO scale. 



5.1.1 Size of the Angular Pixel 

One specific systematic error associated with this method 
is the possible influence that the size of the angular pixel 
has on the determination of the radial BAO scale, since it 
determines which galaxies are considered coUinear in the 
analysis. The effect of different pixel sizes on the radial cor- 
relation function can be seen in Figure [6] The radial correla- 
tion function clearly changes on small scales, but the change 
on the BAO scale is very small. The small scale effect comes 
from the smoothing of the correlation function due to the 
inclusion in the calculation of galaxy pairs which are not ex- 
actly coUinear. For a larger angular pixel, the effect is larger. 
However, the scale where this effect acts is fixed by the an- 
gular pixel size, which is very far from the BAO scale, that 
remains, therefore, unaffected, since the parametrization is 
able to absorb the change in the slope of the function. 

In order to quantify this influence on the determination 
of the BAO scale as a systematic error, we have repeated the 
full analysis for different pixel sizes. The obtained results 
are shown in Figure [7] The radial BAO scale is recovered 
with high precision for any angular pixel size, even for sizes 
as large as 1 degree, which corresponds to a range of two 
orders of magnitude. The associated systematic error can 
be estimated to be S(Azbao) = 0.20%. 



5.1.2 Non-linearities 

Also the error due to the uncertainty in the goodness of 
the description provided by the parametrization for different 
theoretical effects (non-linearities at the scale of the BAO 
peak) has been computed obtaining a global error of 0.10%. 
This was estimated in a conservative way as the difference 
between the Azbao measured using linear and non- linear 
^11 (Az), for the same redshift bins of the analysis. Non- 
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Figure 7. Variation of the radial BAO scale determination as a 
function of the angular pixel size for different redshifts. Results 
are stable, and the maximum variation is always of a few parts 
per mille, very well below 1%, even if the range in pixel sizes 
covers two orders of magnitude. Dashed lines indicate the size 
of the statistical error, including Poisson shot noise and cosmic 
variance, for the used simulation, that covers 1/8 of the sky. 



linearities are computed using the RPT formalism (Crocce 



fc Scoccimarro|2006 1, excluding mode-mode coupling, since 
it only affects small scales, far enough from the BAO scale. 
The contribution of these uncertainties to the systematic 
error can be estimated as S{Azbao) ~ 0.10%. 
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5.1.3 Galaxy Bias 

At the BAO scales it is a good approximation to con- 



sider that the bias is scale independent (Cresswell & Per- 
cival||2069 Crocce et al.||20lT |. It affects the radial correla- 
tion function not only as an overall normalization, but also 
through its effect in the redshift space distortions. Bias can 
influence the determination of the BAO scale only through 
the changes in the goodness of the parametrization of the 
correlation function for different biases. In order to estimate 
the contribution of the galaxy bias to the measurement of 
the radial BAO scale, we have repeated the analysis with 
different values of the bias, to obtain the propagation of the 
uncertainty in the galaxy bias for the selected galaxy pop- 
ulation to the measured value of the radial BAO scale. The 
influence on the peak position is small, and we can estimate 
the associated systematic error as 5{Azbao) = 0.15%. 

Moreover, we have tested the effect of a scale depen- 
dent bias, introducing artificially the effect in the correla- 
tion functions, using an approximated Q-model with the 
determination of parameters of |Cresswell fc Percival| ( [2009[ ). 
The bias variation with A2 in the fitted region of 5||(Aa) 
ranges from 1% to 6%, but the measurement of the BAO 
scale is insensitive to these changes. We estimate the sys- 
tematic error in the presence of a scale dependent bias as 
S{AZBA0) = 0.20%. 
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5.1.4 Starting and End Point of the Fit 

To compute the systematic error associated to the 
parametrization method, we have done some further analy- 
sis on theoretical radial correlation functions with the same 
bin widths and central redshifts as those used in the anal- 
ysis of the MICE simulation. The error associated to the 
method comes from the possible influence in the obtained 
/S.ZBAO of the range of used to perform the fit. To eval- 
uate the error, we have varied this range for the 3 redshift 
bins where we have a significant detection of the BAO scale, 
and performed the fit for each range. 

In the decision of the range to be fitted, we have to 
choose a starting point at angles smaller than the BAO peak, 
where the physics is determined by non-linearities, and an 
end point after the peak, beyond which the effects of cosmic 
variance may be relevant. By varying these two points we 
can study how much the result varies with this decision. 
Results can be seen in Fig. |8] where the obtained Azsao is 
shown for different starting points and end points of the fit, 
for the 3 redshift bins. In all cases, the uncertainty is of the 
order of 0.1%, which we assign as the associated systematic 
error. 



5.1.5 Total Systematic Error 

The different sources of the systematic errors are completely 
independent, and therefore, we can compute the total sys- 
tematic error by summing quadratically these contributions, 
resulting on a value of &sys{^zbao) ~ 0.33%. 

There are some other potential systematic errors, the 
gravitational lensing magnification, which introduces a small 
correlation between redshift bins, or those mainly associated 
to the instrumental effects which could affect the used galaxy 



0.75 <z< 1.05 



1.05<z<1.20 



1.05<z<1.20 



Figure 8. I^zbao evolution as a function of the starting and 
end point of the fitted region. Results are stable, confirming that 
the systematic error is small. Dashed lines indicate the size of the 
statistical error, including Poisson shot noise and cosmic variance, 
for the used simulation, that covers 1/8 of the sky. 
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5.3 



Combination of Radial and Angular BAO 
Scales 



Figure 9. Contours at Icr, 2(t and 3cr C. L. on the plane (wo, Wa) 
(top) and on the plane (Oj\/,uio) (bottom) obtained from the 
analysis of the radial BAO scale. The dot shows the value of the 
parameters for the MICE cosmology. No combination with any 
other cosmological probe is included. The other parameters have 
been fixed to the values of the MICE cosmology. 



We have applied the method of Sanchez et al. (2011 1 to de- 



termine the angular BAO scale for the same simulation. The 
constraints on the {wo,Wa) plane coming from this determi- 
nation of the angular BAO scale can be seen in Figure [TO] 
(top) as thin dashed lines, and correspond to the la, 2a and 
3a C. L. contours. The sensitivity of the angular BAO scale 
is complementary to that of the radial BAO, shown as the 
thin solid lines. When combined, the contours represented by 
the thick solid lines are found. The same constraints for the 
{Qm, Wo) plane are shown in the bottom panel of Figure 10 
These constraints are much stronger than what is usually 
quoted for the BAO standard ruler, which is based on the 
use of the monopole of the 3-D correlation function and does 
not exploit the full sensitivity of the data. It is important to 
remark that these constraints are obtained ONLY with the 
BAO standard ruler, independently of any other cosmologi- 
cal probe, which shows the real power of the standard ruler 
method when the full information is used. 

In order to fairly compare the sensitivity of the present 
method with previous results, we provide also the combined 
result of BAO, both radial and angular, with the distance 
measurementes from CMB using the WMAP7 covariance 



matrix (Komatsu et al. 20111 and assuming the measure- 



ments correspond to the MICE cosmology. The combination 
has been performed following the procedure as detailed in 
Komatsu et al. (20091. The corresponding contours at la, 



2a and 3a C. L. are presented in Figure [TT] both in the 
plane {wo,Wa) (top) and in the plane {flM,wo) (bottom) as 
thick solid lines, and compared with the result using only 
BAO (both radial and angular), which is presented as the 
thin lines. As before, the other parameters have been kept 
fixed. There is an important improvement in the precision of 
the determination of the corresponding parameters in both 
cases, larger in the (HmjWo) plane, showing that a precise 
measurement can be achieved when all the information pro- 
vided by the BAO scale is included in the fit. 



sample. However, these effects are expected to be very small 
and we have neglected them in this analysis. 



5.2 Cosmological Constraints 

The evolution of the measured radial BAO scale, including 
the systematic errors, with redshift is shown in Figure [5] 
The cosmological model of the simulation is the solid line. 
The recovered BAO scale is perfectly compatible with the 
true model, demonstrating that the method works. 

When these measurements are translated into con- 
straints on the cosmological parameters, we obtain the re- 
sults depicted in Figure [9j where the contours for la, 2a 
and 3a C. L. in the {wo,'Wa) plane are shown at the top 
panel and the bottom panel shows the same contours in the 
{Qm,wo) plane. All other parameters have been kept fixed 
to their values for the simulation. The cosmology of the sim- 
ulation is recovered, and the plot shows the sensitivity of the 
radial BAO scale alone, since no other cosmological probe is 
included in these constraints. 



6 CONCLUSIONS 

We have developed a new method to measure the BAO scale 
in the radial two-point correlation function. This method is 
adapted to the observational characteristics of galaxy sur- 
veys, where only the angular position on the sky and the 
redshift are measured for each galaxy. The sound horizon 
scale can be recovered from the non-linear radial correlation 
functions to a very high precision, only limited by the volume 
of the considered survey, since the systematic uncertainties 
associated to the determination of the BAO scale are very 
small, around 0.3%. On the other hand, the method is fully 
cosmology independent, since it relies only on observable 
quantities and, consequently, its results can be analyzed in 
any cosmological model. 

The method has been tested with a mock catalogue 
built upon a large N-body simulation provided by the MICE 
collaboration, in the light cone and including redshift space 
distortions. The true cosmology is recovered within 1-a. An 
evaluation of the main systematic errors has been included 
in this study, and we find that the method is very promis- 
ing and very robust against systematic uncertainties. Note 
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Figure 10. Contours at la, la and Scr C. L. on the plane (uiq, tOa) 
(top) and on the plane (SImi 'i'o) (bottom) from radial BAO (thin 
solid lines), angular BAO (thin dashed lines) and the combination 
of both (thick solid lines) . The dot shows the value of the param- 
eters for the MICE cosmology. No other cosmological probe is 
included in this result, showing the high sensitivity that the BAO 
standard ruler can achieve. The other parameters have been fixed 
to their values in the MICE cosmology. 




Angular + Radial BAO 

-0.6- 

Angular + Radial BAO +WMAP7 

-0.8- 




-1.2 - 



-1.4- 

' ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I 

0.1 0.2 0.3 0.4 



Figure 11. Contours at Icr, 2a and Scr C. L. on the plane (uiq, Wa) 
(top) and on the plane (fijvfi'i'o) (bottom) from BAO (both ra- 
dial and angular), depicted as thin solid lines, and adding also 
the CMB measurements, depicted as thick solid lines. The dot 
shows the value of the parameters for the MICE cosmology. The 
covariance matrix of WMAP7 has been used, while the central 
value of the measurement has been taken at the MICE cosmol- 
ogy. The other parameters have been fixed to their values in the 
MICE cosmology. 



that this analysis over the MICE simulations is done on dark 
matter particles, instead of galaxies. We believe this simplifi- 
cation is not an essential limitation to the method presented 
here, as we have shown that both the modeling and the error 
analysis are quite generic. 

An important conclusion from this study is the achiev- 
able precision on the determination of the cosmological pa- 
rameters that the BAO standard ruler method allows. It 
is much stronger than is usually quoted. The main reason 
is that the current analyses only use the monopole of the 
3-D galaxy correlation function, which does not contain as 
much information as the combination of the angular and ra- 
dial BAO scale. When radial and angular BAO are combined 
the constraints on the cosmological parameters become very 
strong, even without combining them with the CMB. 
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APPENDIX A: THEORETICAL CALCULATION 
OF THE COVARIANCE MATRIX FOR THE 3-D 
CORRELATION FUNCTION 

We have computed a theoretical estimation of the covari- 
ance matrix for the 3-D correlation function, which is then 
applied to the estimation of the covariance for the radial cor- 
relation function. It uses a similar approach to that used in 



Crocce et al. (2011 1 for the angular correlation function. We 



where fik is the cosine of the angle between the line of sight 
and k and Li are the Legendre polynomials. Given a real- 
ization of the density field Sk in a region of volume V, we 
define the estimator for Pi{k) as 



(A3) 



in close analogy with the angular case. 

Assuming that 5k is Gaussianly distributed we can ob- 
tain the covariance for Pi{k): 



{PiAki)PiAk2)) ~ {PiAki)){Pi,{k2)) 



V 47rl/ 



3\ 2 



dQ.kdQ,k,{{5k5k') - {Sk){Sk'))Li^i^J.ki)Li2{fik2)- 



Using Wick's theorem and equation |AT] this covariance can 
be written as 



V 



Y,PiAki)PiAk2)Bi,i^i,i, (A4) 



where 



„ (2Zi-H) (2/2 + 1) f\ , , , , 



which simplifies for the Kaiser approximation to 



47r2 5^(fcl-fc2)„2 



V 



kl 



with 



1-3,Ia 



(A5) 



(A6) 



and the only non-zero coefficients are i'o(/3) = + , 

V2{P) = tP+ 1/3" and Vi{P) = 

The correlation function can also be expanded in Leg- 
endre polynomials. 



(A7) 



have used the following convention for the power spectrum: 



where r is the pair separation and ^ is the cosine of the 
angle between the separation vector and the line of sight. 
The multipoles ^; are directly related to the power spectrum 
multipoles via 



{SkSk') = 5°(fe + k')Pk = ■j^.S^^k'Pk 



(Al) 



(27r)» 

where 5^ and 5^ are the Dirac and Kronecker deltas, and 
V is the volume. In general, the power spectrum can be 
expanded in multipoles. 



Pfe = X^Pi(fc)L,(Mfe) 



(A2) 



6(0 



2^ 



k^ji{kr)Pi{k)dk 



(A8) 



where ji are the spherical Bessel functions. We can, then, 
relate the covariance of the power spectrum with the covari- 
ance of the correlation function. The final result is 



(Ae,,(ri)ACi,(r2)) = 
Ai.iAP) I PkjiAkri)ji,{kr2)k^dk 



(A9) 
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and finally the covariance for the correlation function is 
(Ae(ri,/xi)A^(r2,M2)> = ^{A^iAri)AiiAr2))LiMLiM. 

h,l2 

The covariance for the radial correlation function is obtained 
for the value of /i = 1 which corresponds to the line of sight 
direction. 



